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1 Introduction 



The concept of random dynamical system is a com- 
paratively recent development combining ideas and 
methods from the well developed areas of probabil- 
ity theory and dynamical systems. 

Let us consider a mathematical model of some physi- 

k 

cal process given by the iterates T = Too • • • o7b, k > 
1, of a smooth transformation 7b : M of a man- 
ifold into itself. A realization of the process with 
initial condition xq is modelled by the sequence 
(T£(x ))k>i, the orbit of xq. 

Due to our inaccurate knowledge of the particular 
physical system or due to computational or theo- 
retical limitations (lack of sufficient computational 
power, inefficient algorithms or insufficiently devel- 
oped mathematical or physical theory, for example), 
the mathematical models never correspond exactly to 
the phenomenon they are meant to model. Moreover 
when considering practical systems we cannot avoid 
either external noise or measurement or inaccuracy 
errors, so every realistic mathematical model should 
allow for small errors along orbits not to disturb too 
much the long term behavior. To be able to cope 
with unavoidable uncertainty about the "correct" pa- 



rameter values, observed initial states and even the 
specific mathematical formulation involved, we let 
randomness be embedded within the model to begin 
with. 

We present the most basic classes of models in what 
follows, then define the general concept and present 
some developments and examples of applications. 



2 Dynamics with noise 

To model random perturbations of a transformation 
To we may consider a transition from the image To (x) 
to some point according to a given probability law, 
obtaining a Markov Chain, or, if To depends on a pa- 
rameter p, we may choose p at random at each iter- 
ation, which also can be seen as a Markov Chain but 
whose transitions are strongly correlated. 

2.1 Random noise 

Given T : M Q and a family {p(- \ x) : x € M} 
of probability measures on M such that the sup- 
port of p(- | x) is close to Tq{x), the random or- 
bits are sequences l where each Xk+\ is a ran- 
dom variable with law p{- \ Xk)- This is a Markov 
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Chain with state space M and transition probabil- 
ities {/?(■ | x)} xe M- To extend the concept of in- 
variant measure of a transformation to this setting, 
we say that a probability measure n is stationary if 
fj(A) = J p(A\ x)djx{x) for every measurable (Borel) 
subset A. This can be conveniently translated by say- 
ing that the skew -product measure ^x/onMx M N 
given by 

d(/ux p m )(xq,xi,...,x„,...) 

= d/u(x )p(dxi | x {) ) ■ ■ ■ p(dx n+ i |* n )--- 

is invariant by the shift map S : M x M N O on the 
space of orbits. Hence we may use the Ergodic The- 
orem and get that time averages of every continuous 
observable cp : M — > R, i.e. writing x = (xk)k>o and 

Y n—l 

cp(x) = lim -£(p(x fc ) 



Example 1 (Random jumps). Given £ > and 7b 
M — ► Af, let us define 



>t=o 

J B— 1 



lim - £ <p(n Q (S k (x))) 



>t=o 

exist for x /? N -almost all sequences x, where Jto : 
M x M N — > M is the natural projection on the first 
coordinate. It is well known that stationary measures 
always exist if the transition probabilities p(- \ x) de- 
pend continuously on x. 

A function cp : M — > R is invariant if cp(x) = 
/ (p(z)p(dz | x) for /^-almost every x. We then say 
that /u is ergodic if every invariant function is con- 
stant /^-almost everywhere. Using the Ergodic Theo- 
rem again, if ^ is ergodic, then cp = fq>dp., /^-almost 
everywhere. 

Stationary measures are the building blocks for more 
sophisticated analysis involving e.g. asymptotic so- 
journ times, Lyapunov exponents, decay of correla- 
tions, entropy and/or dimensions, exit/entrance times 
from/to subsets of M, to name just a few frequent 
notions of dynamical and probabilistic/statistical na- 
ture. 



p e (A\x) 



m(AnB(T Q (x),e)) 



m(B(T Q (x),e)) 

where m denotes some choice of Riemannian volume 
form on M. Then p e (- \x) is the normalized volume 
restricted to the e-neighborhood of Tq(x). This de- 
fines a family of transition probabilities allowing the 
points to "jump" from T$(x) to any point in the 8- 
neighborhood of Tq(x) following a uniform distribu- 
tion law. 

2.2 Random maps 

Alternatively we may choose maps T\ , T2, . . . , T\ in- 
dependently at random near To according to a proba- 
bility law v on the space T(M) of maps, whose sup- 
port is close to To in some topology, and consider 
sequences x^ = T^ o ■ ■ ■ o T\ [xo) obtained through ran- 
dom iteration, k > 1 , xq € M. 

This is again a Markov Chain whose transition prob- 
abilities are given for any x E M by 

p(A\x) = v({T eT(M) :T(x) eA}), 

so this model may be reduced to the first one. How- 
ever in the random maps setting we may associate 
to each random orbit a sequence of maps which are 
iterated, enabling us to use robust properties of the 
transformation To (i.e. properties which are known 
to hold for To and for every nearby map T) to derive 
properties of the random orbits. 

Under some regularity conditions on the map x 1— > 
p(A I x) for every Borel subset A, it is possible to 
represent random noise by random maps on suit- 
ably chosen spaces of transformations. In fact the 
transition probability measures obtained in the ran- 
dom maps setting exhibit strong spatial correlation: 
p(- \ x) is close to p(- j y) is x is near y. 
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If we have a parameterized family T : 11 x M — ► M 
of maps we can specify the law v by giving a proba- 
bility 6 on 11 . Then to every sequence T\ , . . . , 7^, . . . 
of maps of the given family we associate a sequence 
GOi , . . . , (Ok, . . . of parameters in 11 since 

r*o...o7i =r fflt o-or ffll m 

for all k > 1, where we write T a (x) = T((0,x). In 
this setting the shift map S becomes a skew-product 
transformation 

S:Mxu N (*,©) i-> (r mi (x),a(ro)), 

to which many of the standard methods of dynami- 
cal systems and ergodic theory can be applied, yield- 
ing stronger results that can be interpreted in random 
terms. 

Example 2 (Parametric noise). Let T : P x M — ► M 
be a smooth map where P,M are finite dimensional 
Riemannian manifolds. We fix po € P, denote by 
m some choice of Riemannian volume form on P, 
set T w {x) = T(w,x) and for every £ > write e = 
(m(B(po,e))~ l ■ (m | B(po,e)), the normalized re- 
striction of m to the e-neighborhood of po. Then 
(T w )weP together with e defines a random pertur- 
bation of T po , for every small enough £ > 0. 

Example 3 (Global additive perturbations). Let M 
be a homogeneous space, i.e., a compact connected 
Lie Group admitting an invariant Riemannian metric. 
Fixing a neighborhood 11 of the identity e € M we 
can define a map T : 11 xM —>M,(u,x) i— > L u (Tq{x)), 
where L u (x) = u ■ x is the left translation associated 
to w £ M. The invariance of the metric means that 
left (an also right) translations are isometries, hence 
fixing u G 11 and taking any (x, v) € TM we get 

\\DT u (x)-v\\ = \\DL u (T Q (x))(DT (x)-v)\\ 
= \\DT (x)-v\\. 

In the particular case of M = T d , the J-dimensional 
torus, we have T u (x) = Tq(x) + u and this simplest 



case suggests the name additive random perturba- 
tions for random perturbations defined using families 
of maps of this type. 

For the probability measure on 11 we may take 
9 8 any probability measure supported in the £- 
neighborhood of e and absolutely continuous with 
respect to the Riemannian metric on M, for any £ > 
small enough. 

Example 4 (Local additive perturbations). If M = 
W 1 and Uq is a bounded open subset of M 
strictly invariant under a diffeomorphism 7b, i.e., 
closure (To(Uo)) C Uq, then we can define an isomet- 
ric random perturbation setting 

• V = T (Uq) (so that closure (V) = 
closure {T (U )) C U ); 

• G ~ M. d the group of translations of R d ; 

• V a small enough neighborhood of in G. 

Then for v G V and iGVwe set T v (x) = x + v, with 
the standard notation for vector addition, and clearly 
T v is an isometry. For 6 e we may take any probability 
measure on the £-neighborhood of 0, supported in <V 
and absolutely continuous with respect to the volume 
in M. d , for every small enough £ > 0. 

2.3 Random perturbations of flows 

In the continuous time case the basic model to 
start with is an ordinary differential equation dX t = 
f(t,X t )dt, where / : [0, +°o) -► X (Af) and X (M) is 
the family of vector fields in M. We embed random- 
ness in the differential equation basically through 
diffusion, the perturbation is given by white noise or 
Brownian motion "added" to the ordinary solution. 

In this setting, assuming for simplicity that M = W 1 , 
the random orbits are solutions of stochastic differ- 
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ential equations 



2.4 The abstract framework 



dXt = fit , X, ) dt + £ • c(f , X t ) dW t , < t < T, X = Z, 

where Z is a random variable, £, T > and both / : 
[0, T] x R" -» R" and a : [0, r] x M" -> £ (R* R") are 
measurable functions. We have written x(R*,R") 
for the space of linear maps R k — > R" and for the 
w/u'fe «o«e process on R*. The solution of this equa- 
tion is a stochastic process 

I:lx(l^M, (t,(0)^X t ((0), 

for some (abstract) probability space £2, given by 

X ; = Z+ f T f(s,X s )ds+ f J e-a{s,X s )dW s , 
Jo Jo 

where the last term is a stochastic integral in the 
sense of Ito. Under reasonable conditions on / and 
a, there exists a unique solution with continuous 
paths, i.e. 

[0,+°o) 3t^X t ((Q) 

is continuous for almost all co G £2 (in general these 
paths are nowhere differentiable). 

Setting Z = 8^ the probability measure concentrated 
on the point xo, the initial point of the path is xq with 
probability 1. We write X t ((o)xo for paths of this 
type. Hence x i— > X t (oo)x defines a map X t (co) : M O 
which can be shown to be a homeomorphism and 
even a diffeomorphisms under suitable conditions on 
/ and a. These maps satisfy a cocycle property 

Xo(co) = Mm (identity map of M), 
X t+S (a) = X t (e(s){a)) o^(co), 

for s,t > and co G Q, for a family of measure pre- 
serving transformations Q(s) : (il,P) O on a suitably 
chosen probability space (Q,P). This enables us to 
write the solution of this kind of equations also as a 
skew-product. 



The illustrative particular cases presented can all be 
written in skew-product form as follows. 

Let (£2,P) be a given probability space, which will 
be the model for the noise, and T be time, which 
usually means Z + ,Z (discrete, resp. invertible sys- 
tem) or R + ,R (continuous, resp. invertible system). 
A random dynamical system is a skew-product 

5,:flxMO, (co,*) i-> (6(f)(eo),<p(f,eo)(jc)), 

for all t G T, where 6 : T x ^ is a family of 
measure preserving maps 8(f) : (Q,P) O and (p : 
TxdxM^Misa family of maps cp(f , co) : M O 
satisfying the cocycle property: for s,t G T, co G £2 

cp(0,co) = ld M , 
(p(f + s,co) = q>(f, 8 ($)(©)) o(p(j,co). 

In this general setting an invariant measure for the 
random dynamical system is any probability measure 
p. on Q. x M which is 5 r invariant for all f G T and 
whose marginal is P, i.e. fi(S,-\U)) =fi(U) and 
h(tIq (U)) = ¥(U) for every measurable U C £1 x 
M, respectively, with jiq : Q. x M — > £1 the natural 
projection. 

Example 5. In the setting of the previous examples of 
random perturbations of maps, the product measure 
Tl = Px i t/on£lxM, with £2 = U N , F = 8^ and,u any 
stationary measure, is clearly invariant. However not 
all invariant measures are product measures of this 
type. 

Naturally an invariant measure is ergodic if every S t - 
invariant function is /v-almost everywhere constant, 
i.e. if \\f : £2 x M — > R satisfies \\tos t = Y /^-almost 
everywhere for every f G T, then \|/ is /^-almost ev- 
erywhere constant. 



4 



3 Applications 

We avoid presenting well established applications of 
both probability or stochastic differential equations 
(solution of boundary value problems, optimal stop- 
ping, stochastic control etc) and dynamical systems 
(all sort of models of physical, economic or bio- 
logical phenomena, solutions of differential equa- 
tions, control systems etc), focusing instead on top- 
ics where the subject sheds new light on these areas. 

3.1 Products of random matrices and the 
Multiplicative Ergodic Theorem 

The following celebrated result on products of ran- 
dom matrices has far-reaching applications on dy- 
namical systems theory. 

Let (X„) n >o be a sequence of independent and iden- 
tically distributed random variables on the probabil- 
ity space (£1,P) with values in L (M. k ,R k ) such that 
£"(log + H-Xill) < +00, where log + x = max{0,logx} 
and || • || is a given norm on £(M*,]R*). Writing 
(Pn(co) = X ra (co) o • • • 0X1 (go) for all n > 1 and GO <E £1 
we obtain a cocycle. If we set 



B 



a' 



1 



lim -log||(p„(C0)y|| 



= {(co,y) eflxK': 

exists and is finite or is — 00} ; and 
= {G0GQ: (co,y) for ally £R k }, 



then Q.' contains a subset Q." of full probability and 
there exist random variables (which might take the 
value — 00) X,i > X2 > ■ • • > 'kk with the following 
properties. 

(1) Let / = {k + 1 = ?'i > ?2 > ■ ■ • > = 1 } be any 
(Z + 1) -tuple of integers and then we define 

£2/ = {go € Q." : X,,(go) = Xj((£t),ih > i,j > h+\, and 
h h (go) > h h+l (co) for all 1< h < 1} 



the set of elements where the sequence X, jumps ex- 
actly at the indexes in /. Then for GO G £2/, 1 <h <l 



^((d) = {y e : 



lim -log 1 1 q>„ (go) 



is a vector subspace with dimension i^-i 
(2) Setting E /; * +1 (go) = {0}, then 

1 



lim -log 1 1 q>„ (00) 



for every y € £/,/,(a>) \ Z/^+i (00). 

(3) For all go <G Q." there exists the matrix 



A(a>) 



lim 



(cp„(G0))*(p„(G0) 



l/2n 



whose eigenvalues form the set {e^ :i= 1, . . . ,k}. 

The values of X; are the random Lyapunov charac- 
teristics and the corresponding subspaces are analo- 
gous to random eigenspaces. If the sequence (X„)„>o 
is ergodic, then the Lyapunov characteristics become 
non-random constants, but the Lyapunov subspaces 
are still random. 

We can easily deduce the Multiplicative Ergodic 
Theorem for measure preserving differentiable maps 
(To, ft) on manifolds M from this result. We as- 
sume for simplicity that M C R k and set p(A\x) = 
h To i x \(A) = 1 if Tq(x) € A and otherwise. Then the 
measure /ux p N onM x M N is a-invariant (as defined 
in Section 2) and we have that Jio 00 = 7o o7t , where 
Jlo : M N — > M is the projection on the first coordinate, 
and also (7to)*(u x p N ) = fj. Then setting for n > 1 



X 



M 

x 



L 



DT Q (x) 



and X n = X o jt o o n 



we obtain a stationary sequence to which we can ap- 
ply the previous result, obtaining the existence of 
Lyapunov exponents and of Lyapunov subspaces on 
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a full measure subset for any C measure preserving 
dynamical system. 

By a standard extension of the previous setup we ob- 
tain a random version of the multiplicative ergodic 
theorem. We take a family of skew-product maps 
St '■ Q x M O as in Subsection 12.41 with an invariant 
probability measure p and such that <p(f,G)) : M 
is (for simplicity) a local diffeomorphism. We then 
consider the stationary family 

X t : Q. L{TM) 

CO h-> Dq>(f,a>) : TM Q ' 

where Dcp(7,co) is the tangent map to cp(?, to). This is 
a cocycle since for all t,s G T, CO € 12 we have 

X(j+f,co) =X(j,e(f)ffl) oX(/,co). 

If we assume that 

sup sup (log + ||Dcp(?,co)(x)||) £L l (Q,F), 

0<t<\ x€M 

where || • || denotes the norm on the corresponding 
space of linear maps given by the induced norm 
(from the Riemannian metric) on the appropriate tan- 
gent spaces, then we obtain a sequence of random 
variables (which might take the value — oo) X\ > X2 > 
• • • > Xk, with k being the dimension of M, such that 

lim -log\\X t (a>,x)y\\ =Xi((0,x) 
t— >+<*> t 

for every y € £,a),x) = Z,(co,x) \E/+i (co,x) and i = 
l,...,k+l where (£,(o),x)); is a sequence of vector 
subspaces in T X M as before, measurable with respect 
to (co,x). In this setting the subspaces Ef((0,x) and 
the Lyapunov exponents are invariant, i.e. for all t G 
T and p almost every ((0,x) GflxMwe have 

Xi(S t ((Q,x)) = Xi((H,x) and £' ! -(5r ((0,x)) =£ , ( -((0,x). 

The dependence of Lyapunov exponents on the map 
Tq has been a fruitful and central research program 



in dynamical systems for decades extending to the 
present day. The random multiplicative ergodic the- 
orem sets the stage for the study of the stability of 
Lyapunov exponents under random perturbations. 

3.2 Stochastic stability of physical measures 

The development of the theory of dynamical sys- 
tems has shown that models involving expressions 
as simple as quadratic polynomials (as the logistic 
family or Henon attractor), or autonomous ordinary 
differential equations with a hyperbolic singularity 
of saddle-type, as the Lorenz flow, exhibit sensitive 
dependence on initial conditions, a common feature 
of chaotic dynamics: small initial differences are 
rapidly augmented as time passes, causing two tra- 
jectories originally coming from practically indistin- 
guishable points to behave in a completely different 
manner after a short while. Long term predictions 
based on such models are unfeasible since it is not 
possible to both specify initial conditions with arbi- 
trary accuracy and numerically calculate with arbi- 
trary precision. 

Physical measures 

Inspired by an analogous situation of unpre- 
dictability faced in the field of Statistical Mechan- 
ics/Thermodynamics, researchers focused on the 
statistics of the data provided by the time averages of 
some observable (a continuous function on the man- 
ifold) of the system. Time averages are guaranteed 
to exist for a positive volume subset of initial states 
(also called an observable subset) on the mathemat- 
ical model if the transformation, or the flow associ- 
ated to the ordinary differential equation, admits a 
smooth invariant measure (a density) or a physical 
measure. 
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Indeed, if fuo is an ergodic invariant measure for the 
transformation To, then the Ergodic Theorem ensures 
that for every ^-integrable function cp : M — > R and 
for /^-almost every point x in the manifold M the 
time average cp(x) = lim n ^ +00 n~ 1 cp(r o 7 (x)) ex- 
ists and equals the space average / cpii/vo- A phys- 
ical measure p. is an invariant probability measure 
for which it is required that time averages of every 
continuous function (p exist for a positive Lebesgue 
measure (volume) subset of the space and be equal 
to the space average ju(cp). 

We note that if /j is a density, that is, is abso- 
lutely continuous with respect to the volume mea- 
sure, then the Ergodic Theorem ensures that p. is 
physical. However not every physical measure is ab- 
solutely continuous. To see why in a simple example 
we just have to consider a singularity p of a vector 
field which is an attracting fixed point (a sink), then 
the Dirac mass 5 P concentrated on p is a physical 
probability measure, since every orbit in the basin of 
attraction of p will have asymptotic time averages for 
any continuous observable cp given by cp(/?) = 8 p ((p). 

Physical measures need not be unique or even ex- 
ist in general, but when they do exist it is desirable 
that the set of points whose asymptotic time aver- 
ages are described by physical measures (such set 
is called the basin of the physical measures) be of 
full Lebesgue measure — only an exceptional set of 
points with zero volume would not have a well de- 
fined asymptotic behavior. This is yet far from being 
proved for most dynamical systems, in spite of much 
recent progress in this direction. 

There are robust examples of systems admitting sev- 
eral physical measures whose basins together are 
of full Lebesgue measure, where robust means that 
there are whole open sets of maps of a manifold in 
the C 2 topology exhibiting these features. For typical 
parameterized families of one-dimensional unimodal 



maps (maps of the circle or of the interval with a 
unique critical point) it is known that the above sce- 
nario holds true for Lebesgue almost every parame- 
ter. It is known that there are systems admitting no 
physical measure, but the only known cases are not 
robust, i.e. there are systems arbitrarily close which 
admit physical measures. 

It is hoped that conclusions drawn from models ad- 
mitting physical measures to be effectively observ- 
able in the physical processes being modelled. In 
order to lend more weight to this expectation re- 
searchers demand stability properties from such in- 
variant measures. 

Stochastic stability 

There are two main issues when we are given a 
mathematical model, both theoretical but with prac- 
tical consequences. The first one is to describe the 
asymptotic behavior of most orbits, that is, to un- 
derstand where do orbits go when time tends to in- 
finity. The second and equally important one is to 
ascertain whether the asymptotic behavior is stable 
under small changes of the system, i.e. whether the 
limiting behavior is still essentially the same after 
small changes to the evolution law. In fact since 
models are always simplifications of the real system 
(we cannot ever take into account the whole state of 
the universe in any model), the lack of stability con- 
siderably weakens the conclusions drawn from such 
models, because some properties might be specific to 
it and not in any way resemblant of the real system. 

Random dynamical systems come into play in this 
setting when we need to check whether a given 
model is stable under small random changes to the 
evolution law. 

In more precise terms, we suppose that we are given 
a dynamical system (a transformation or a flow) ad- 
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mitting a physical measure /Jo, and we take any 
random dynamical system obtained from this one 
through the introduction of small random perturba- 
tions on the dynamics, as in Examples [H S or in 
Subsection 12.31 with the noise level £ > close to 
zero. 

In this setting if, for any choice fj e of invariant mea- 
sure for the random dynamical system for all £ > 
small enough, the set of accumulation points of the 
family (ju e ) e> o, when £ tends to — also known as 
zero noise limits — is formed by physical measures 
or, more generally, by convex linear combinations of 
physical measures, then the original unperturbed dy- 
namical system is stochastically stable. 

This intuitively means that the asymptotic behavior 
measured through time averages of continuous ob- 
servables for the random system is close to the be- 
havior of the unperturbed system. 

Recent progress in one-dimensional dynamics has 
shown that, for typical families (/f)/e(o,i) of maps of 
the circle or of the interval having a unique critical 
point, a full Lebesgue measure subset T of the set of 
parameters is such that, for t G T, the dynamics of f t 
admits a unique stochastically stable (under additive 
noise type random perturbations) physical measure 
/ii t whose basin has full measure in the ambient space 
(either the circle or the interval). Therefore models 
involving one-dimensional unimodal maps typically 
are stochastically stable. 

In many settings (e.g. low dimensional dynamical 
systems) Lyapunov exponents can be given by time 
averages of continuous functions — for example the 
time average of log ||ZD7b|| gives the biggest expo- 
nent. In this case stochastic stability directly implies 
stability of the Lyapunov exponents under small ran- 
dom perturbations of the dynamics. 

Example 6 (Stochastically stable examples). Let 7b : 
S 1 O be a map such that X, the Lebesgue (length) 



measure on the circle, is To-invariant and ergodic. 
Then X is physical. 

We consider the parameterized family T t : S 1 x S 1 — > 
S 1 ,^,*) i— >x + t and a family of probability mea- 
sures 8 e = (k(-e,e)y l ■ (X | (-£,£)) given by the 
normalized restriction of X to the £-neighborhood of 
0, where we regard S 1 as the Lie group R/Z and use 
additive notation for the group operation. Since X is 
T, -invariant for every t G § , X is also an invariant 
measure for the measure preserving random system 

for every £ > 0, where CI = Hence (7b, A,) is 

stochastically stable under additive noise perturba- 
tions. 

Concrete examples can be irrational rotations, 
Tq(x) = x + a with a G R \ Q, or expanding maps of 
the circle, Tq{x) = b-x for some b G N, n > 2. Anal- 
ogous examples exist in higher dimensional tori. 

Example 1 (Stochastic stability depends on the type 
of noise). In spite of the straightforward way to ob- 
tain stochastic stability in Example|6l for e.g. an ex- 
panding circle map Tq(x) =2-x, we can choose a 
continuous family of probability measures 6 e such 
that the same map 7b is not stochastically stable. 

It is well known that X is the unique absolutely con- 
tinuous invariant measure for 7b and also the unique 
physical measure. Given £ > small let us define 
transition probability measures as follows 

/ | ) = M H>e(z)-£,4>e(z)+£] 
Pe[ ' lz) HMz)-eMz) + e])' 

where (j) E | (-£,£) = 0, (j) E | [S 1 \ (-2e,2e)] = T and 
over (— 2e, — e] U [e,2e) we define (j) E by interpolation 
in order that it be smooth. 

In this setting every random orbit starting at (—£,£) 
never leaves this neighborhood in the future. More- 
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over it is easy to see that every random orbit eventu- 
ally enters (—£,£). Hence every invariant probabil- 
ity measure /u £ for this Markov Chain model is sup- 
ported in [—£,£]. Thus letting £->0we see that the 
only zero-noise limit is 80 the Dirac mass concen- 
trated at 0, which is not a physical measure for 7q. 

This construction can be done in a random maps set- 
ting, but only in the C° topology — it is not possible 
to realize this Markov Chain by random maps that 
are C 1 close to 7b for £ near 0. 



3.3 Characterization of measures satisfying 
the Entropy Formula 

A lot of work has been put in recent years in extend- 
ing important results from dynamical systems to the 
random setting. Among many examples we mention 
the local conjugacy between the dynamics near a hy- 
perbolic fixed point and the action of the derivative 
of the map on the tangent space, the stable/unstable 
manifold theorems for hyperbolic invariant sets and 
the notions and properties of metric and topological 
entropy, dimensions and equilibrium states for po- 
tentials on random (or fuzzy) sets. 

The characterization of measures satisfying the En- 
tropy Formula is one important result whose exten- 
sion to the setting of iteration of independent and 
identically distributed random maps has recently had 
interesting new consequences back into non-random 
dynamical systems. 



Metric entropy for random perturbations 

Given a probability measure /u and a partition t, of 
M, except perhaps for a subset of ^-null measure, the 



entropy of n with respect to ^ is defined to be 

= -£/!(*) lOg/l(Z0 

Ret, 

where we convention OlogO = 0. Given another fi- 
nite partition C, we write t, V C, to indicate the partition 
obtained through intersection of every element of ^ 
with every element of C,, and analogously for any fi- 
nite number of partitions. If /u is also a stationary 
measure for a random maps model (as in Subsec- 
tion 12.2b . then for any finite measurable partition 
ofM, 

V© = inf \ J H/\/(T^-\^)dp N (^ 

i'=0 

is finite and is called the entropy of the random dy- 
namical system with respect to ^ and to p.. 

We define h u = sup^ h^) as the metric entropy of 
the random dynamical system, where the supremo 
is taken over all ^-measurable partitions. An impor- 
tant point here is the following notion: setting A the 
Borel G-algebra of M, we say that a finite partition £, 
of M is a random generating partition for A if 

i=0 

(except ^-null sets) for /7 N -almost all CO € £1 = 
Then a classical result from Ergodic Theory ensures 
that we can calculate the entropy using only a ran- 
dom generating partition t,: we have h^ = fyu(^). 

The Entropy Formula 

There exists a general relation ensuring that the en- 
tropy of a measure preserving differentiable transfor- 
mation (To,/u) on a compact Riemannian manifold is 
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bounded from above by the sum of the positive Lya- 
punov exponents of Tq 

h,(T Q )< [ £ hix)dn{x). 

J Xi(x)>0 

The equality (Entropy Formula) was first shown 
to hold for diffeomorphisms preserving a measure 
equivalent to the Riemannian volume, and then the 
measures satisfying the Entropy Formula were char- 
acterized: for C 2 diffeomorphisms the equality holds 
if, and only if, the disintegration of fu along the un- 
stable manifolds is formed by measures absolutely 
continuous with respect to the Riemannian volume 
restricted to those submanifolds . The unstable man- 
ifolds are the submanifolds of M everywhere tangent 
to the Lyapunov subspaces corresponding to all pos- 
itive Lyapunov exponents, the analogous to "inte- 
grating the distribution of Lyapunov subspaces cor- 
responding to positive exponents" — this particular 
point is a main subject of smooth ergodic theory for 
non-uniformly hyperbolic dynamics. 

Both the inequality and the characterization of sta- 
tionary measures satisfying the Entropy Formula 
were extended to random iterations of independent 
and identically distributed C 2 maps (non-injective 
and admitting critical points), and the inequality 
reads 

h f- L Xi(x,(o)d^(x)dp N ((o). 

J J %i{x,<a)>0 

where the functions X, are the random variables pro- 
vided by the Random Multiplicative Ergodic Theo- 
rem. 

3.4 Construction of physical measures as 
zero-noise limits 

The characterization of measures which satisfy the 
Entropy Formula enables us to construct physical 



measures as zero-noise limits of random invariant 
measures in some settings, outlined in what follows, 
obtaining in the process that the physical measures 
so constructed are also stochastically stable. 

The physical measures obtained in this manner ar- 
guably are natural measures for the system, since 
they are both stable under (certain types of) random 
perturbations and describe the asymptotic behavior 
of the system for a positive volume subset of ini- 
tial conditions. This is a significant contribution to 
the state-of-the-art of present knowledge on Dynam- 
ics from the perspective of Random Dynamical Sys- 
tems. 

Hyperbolic measures and the Entropy Formula 

The main idea is that an ergodic invariant measure p 
for a diffeomorphism To which satisfies the Entropy 
Formula and whose Lyapunov exponents are every- 
where non-zero (known as hyperbolic measure) nec- 
essarily is a physical measure for Tq. This follows 
from standard arguments of smooth non-uniformly 
hyperbolic ergodic theory. 

Indeed /u satisfies the Entropy Formula if, and only if, 
(i disintegrates into densities along the unstable sub- 
manifolds of To. The unstable manifolds W u (x) are 
tangent to the subspace corresponding to every pos- 
itive Lyapunov exponent at /^-almost every point x, 
they are an invariant family, i.e. Tq(W u (x)) = W"(x) 
for /^-almost every x, and distances on them are uni- 
formly contracted under iteration by T ~ l . 

If we know that the exponents along the complemen- 
tary directions are non-zero, then they must be nega- 
tive and smooth ergodic theory ensures that there ex- 
ist stable manifolds, which are submanifolds W s (x) 
of M everywhere tangent to the subspace of negative 
Lyapunov exponents at /^-almost every point x, form 
a To-invariant family (To(W s (x)) = W s (x), /^-almost 
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everywhere), and distances on them are uniformly 
contracted under iteration by Tq. 

We still need to understand that time averages are 
constant along both stable and unstable manifolds, 
and that the families of stable and unstable manifolds 
are absolutely continuous, in order to realize how an 
hyperbolic measure is a physical measure. 

Given y € W s (x) the time averages of x and y co- 
incide for continuous observables simply because 
dist (r ( j 1 (jc), T n (y)) — > when n — > +00. For unsta- 
ble manifolds the same holds when considering time 
averages for T Q ~ 1 . Since forward and backward time 
averages are equal ^-almost everywhere, we see that 
the set of points having asymptotic time averages 
given by /j has positive Lebesgue measure if the fol- 
lowing set 



B = \J{W s (y) :y ew"(x)nsupp(^)} 



has positive volume in M, for some x whose time 
averages are well defined. 

Now, stable and unstable manifolds are transverse 
everywhere where they are defined, but they are only 
defined ^-almost everywhere and depend measur- 
ably on the base point, so we cannot use transver- 
sality arguments from differential topology, in spite 
of W u (x) Hsupp (ji) having positive volume in W u (x) 
by the existence of a smooth disintegration of /u along 
the unstable manifolds. However it is known for 
smooth (C 2 ) transformations that the families of sta- 
ble and unstable manifolds are absolutely continu- 
ous, meaning that projections along leaves preserve 
sets of zero volume. This is precisely what is needed 
for measure-theoretic arguments to show that B has 
positive volume. 



Zero-noise limits satisfying the Entropy Formula 

Using the extension of the characterization of mea- 
sures satisfying the Entropy Formula for the random 
maps setting, we can build random dynamical sys- 
tems, which are small random perturbations of a map 
To, having invariant measures /u e satisfying the En- 
tropy Formula for all sufficiently small £ > 0. Indeed 
it is enough to construct small random perturbations 
of Tq having absolutely continuous invariant proba- 
bility measures /u e for all small enough £ > 0. 

In order to obtain such random dynamical systems 
we choose families of maps T : u x M —> M and 
of probability measures (0 e ) E >o as in Examples [3] 
and |4] where we assume that o G 11 so that To be- 
longs to the family. Letting T x (u) = T(u,x) for all 
(u,x) G U x M, we then have that r v (8 E ) is abso- 
lutely continuous. This means that sets of perturba- 
tions of positive 6 e -measure send points of M onto 
positive volume subsets of M. This kind of perturba- 
tion can be constructed for every continuous map of 
any manifold. 

In this setting we have that any invariant probabil- 
ity measure for the associated skew-product map 
S : £2 x M O of the form Qf x /v e is such that /j £ is 
absolutely continuous with respect to volume on M. 
Then the Entropy Formula holds 




Having this and knowing the characterization of 
measures satisfying the Entropy Formula, it is nat- 
ural to look for conditions under which we can guar- 
antee that the above inequality extends to any zero- 
noise limit /no of jt/g when £ — > 0. In that case /no sat- 
isfies the Entropy Formula for To. 

If in addition we are able to show that p.o is a hyper- 
bolic measure, then we obtain a physical measure for 
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To which is stochastically stable by construction. 

These ideas can be carried out completely for hy- 
perbolic diffeomorphisms, i.e. maps admitting an 
continuous invariant splitting of the tangent space 
into two sub-bundles E © F defined everywhere with 
bounded angles, whose Lyapunov exponents are 
negative along E and positive along F. Recently 
maps satisfying weaker conditions where shown to 
admit stochastically stable physical measures fol- 
lowing the same ideas. 

These ideas also have applications to the construc- 
tion and stochastic stability of physical measure for 
strange attractors and for all mathematical models 
involving ordinary differential equations or iterations 
of maps. 
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